Production of J/^-meson pairs and 4c-tetraquark at LHC 
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Theoretical predictions for pp — )• 2J/^ + X cross section at ^/s = 7 TeV with 
different kinematical restrictions are presented. Results are compared with the the 
first LHCb data available. Special attention is payed to possible signal from novel 
particles — tetraquarks build from two valence c-quarks and two valence c-quarks. 
According to our estimates it is quite possible to observe at least one of these states 
(tensor tetraquark) experimentally under LHCb conditions. 
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I. INTRODUCTION 

Simultaneous production of two J/?/'- mesons was first observed in 1982 by the NA3 collab- 
oration in multi-muon events in pion-platinum interactions at 150 and 280 GeV ll and later 

n 

at 400 GeV in proton-platinum collisions [2]. Cross-section of this reaction is suppressed 
compared to the single J / ip production by more than three orders of magnitude, because of 
both higher order of perturbative QCD O (aj) and kinematical reasons due to production 
of bigger mass. At NA3 energies main contribution to the cross section arose from quark- 
antiquark annihilation channel On the contrary, it is the gg — ?■ J/tpJ/ip process studied 
in ^ which dominates at the LHC. Both contributions lead to the cross-section decreasing 
as the power of the invariant mass of the J/ip pair (a ~ for the gg-annihilation and 
a ~ 1/s^ for the gluonic production). It is worth mentioning that such decrease is absent 
in higher orders of perturbation theory and the cross-section of the gg — 2J/iIj + X process 
becomes constant at large invariant masses ^] . In current work we will be mainly interested 
in the low invariant mass region as it is most accessible in the LHCb experiment, which 
specializes in heavy quarks studies. 
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In the leading order in quarkonia pair production process obeys selection rules similar 
to those in quarkonia decays. Two initial gluons in a color-singlet state are C-even. That 
is why production of J /ip--, f]c- or Xc-nieson pairs is allowed while combined production of 2 
particles having different C-parity (such as J/i' + rjc and J/ip + Xc) is prohibited. 

First comparison of the NA3 data with QCD predictions has shown 2 times shortage of 
the prediction jsj. Taking additional contribution from the P-wave state Xc or ^/''-meson 
decay to J/ip into account could improve the agreement. Selection rules mentioned above 
put restrictions on the feed down from the higher states. For instance feed down from the 
J/i^ + Xc production is expected to be less than from the Xc + Xc one as production of the 
former final state is prohibited in leading order. Meanwhile feed down from the J/tp + ip' 
channel is expected to be quite large. In principle, other mechanisms of J/tp-paii production 
are possible at the same order in as- First of all it is contribution of the color-octet states 
which plays significant role in the single J/ ip production in the high-pr region. Cross section 
of the color-octet states production falls more slowly as transverse momenta grows than that 
of the color-singlet one. However at small transverse momenta and small invariant masses 
of the J /if) pair smallness of the color-octet matrix element in the wave function of the 
J / ip meson compared to the color-singlet one obviously leads to the negligible contribution. 
As shown in [6|, [D] octet contribution becomes significant only at transverse momentum 
Pt ^ 15 GeV (which corresponds to big invariant mass of the J/ip pair). Another possibility 
to enlarge theoretical prediction is connected with the rearrangement of c-quarks in the 
final state, when instead of two colorless cc-pairs two diquarks, cc and cc, are produced. 
After interaction in the final state they can form a pair of double-heavy baryons above the 
threshold of their production or a pair of J/ip-mesons below it. 

The region of small invariant masses of the J/tp-pair is most interesting due to the op- 
portunity for two diquarks [ccjg^ + [cc]^^ to form a bound state — tetraquark, which decays 
to a J/ip-paii. Attraction between 3c- and 3c-states does not exclude such a_possibility, es- 
pecially since similar exotic states like y(3940) resonance decaying to J/ipu 



resonance decaying to J/ip4> 
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and r(4140) 



, ll| have recently been observed. 
In addition, pair production of J/-?/'- mesons has been earlier discussed as a possible way 
to observe C-even states of bottmonia. Production of scalar and tensor x^-mesons was 



considered in 
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13| . production of an r/^-meson — in 



The second section of our article is devoted to the non- resonant production of J/ip- 
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FIG. 1. Feynman diagrams contributing gg — ;> 2J/ip process. 



meson pairs in the gluon-gluon interaction. In the third section cross section of this process 
at LHC at 7 TeV energy is calculated taking different experimental restrictions into account. 
Special attention is payed to production in LHCb conditions as there is first experimental 



data available 



15| . Fourth section is devoted to calculation of the 4c-tetraquark mass and 



estimation of the cross section of its production at LHC. 



II. DOUBLE J/V-MESON PRODUCTION IN GLUON-GLUON INTERACTION 

In the leading order of perturbative QCD there are 31 Feynman diagrams describing color- 
singlet charmonium pairs production in gluonic reaction (see figH]). We are not considering 
contribution of the quark- ant iquark interaction, which is negligibly small at LHC energies. 
Hadronization of a cc-pair into a final J/-?/;- meson is accounted for by the wave function of 
this particle at origin: 

V'J/^WL=o = 0-21GeV3^2^ (1) 

It is the only nonperturbative parameter in the matrix element of gg — ?■ 2J/iIj process. Its 
value is extracted from the leptonic width of J/ip-meson neglecting QCD corrections as we 
do not take these corrections into account in our matrix element. Color-octet contributions 
in the kinematical region considered in our article can be neglected 6|, l7|. 

In the following we will use two methods to calculate the cross section. First approach 
involves standard procedure of analytical calculation of the amplitude, analytical squaring 



of this amplitude (for example, using FeynCalc package [l6|) and subsequent integration 
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TABLE I. Cross section of gg — >■ 2J/-0 processes for different combinations of final particle's polarizations 



Vi, GeV 


a{LL), pb 


aiLT), pb 


a(rT), pb 


6.2 


1.59 


4.66 


5.75 


6.5 


1.75 


5.04 


6.69 


7.0 


1.47 


3.99 


5.96 


8.0 


0.86 


2.03 


3.75 


9.0 


0.50 


1.09 


2.26 


10.0 


0.30 


0.59 


1.39 


12.0 


0.11 


0.196 


0.57 


14.0 


0.046 


0.074 


0.26 


16.0 


0.022 


0.031 


0.136 


18.0 


0.011 


0.014 


0.071 


20.0 


0.006 


0.0071 


0.043 



over the phase space. Second method which we use is based on numerical calculation o " 
the amplitude at each point of the phase space, followed by squaring (see papers [l^, [l8| 
for more details). The latter method is associated with a smaller number of analytical 
calculations and makes it easier to consider different quantum numbers of the final particles 
and their polarizations. Comparison of results obtained by these two methods also help to 
check correctness of the calculations. Both methods give values that agree with each other 
and coincide within errors of calculations with results presented in {4, 19 1. 



In figj2]the cross section of double JZ-^-meson production in gluonic interaction versus 
invariant mass of the pair is shown for different combinations of polarizations. Table |T] shows 
cross section of the J/ip pair production for various polarizations of each of J/'i/'-particles. It 
can be seen that near the threshold cross sections of transversely and longitudinally polarized 
J/ t/'- meson production are comparable, while at higher pt (and hence large invariant mass 
of the gluon pair), J/-?/'- mesons are mainly transversely polarized (see also As an 

example, in figj3] we show angular distribution of polarized J/ ^/'-mesons, produced in gg 
2J/%jj reaction at = 10 GeV. One can easily see, that angular distributions for different 
combinations of polarizations differ drastically. 
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FIG. 2. Cross section of double J/^-meson gluonic production as a function of their invariant mass: solid 
curve — total cross section; dashed curve — both mesons are transversely polarized; dotted curve — one me- 
son is transversely polarized, another — longitudinally; dash-dotted curve — both mesons are longitudinally 
polarized. 



Partonic cross-section summed over polarizations as a function of s is equal 

5127r3^(0)V 
ois) = / X 



X I s(5s^ + 874m^s^ + 13368m^s^ + 36594m'*s + 90m^) \l 1 - ^+ 



1 + a/I 



+ 30m^(6m^ + 70m^s - 1025m^s^ - 134s^) log 



4m2 



1- a/1- 



(2) 



Near the threshold there is a typical square- root dependence (a ~ y's — 4mj^^), which 
corresponds to the production of a pair in S'-wave state. At large invariant masses expression 
(E]) tends to ~ s"^ asymptotics. 

As noted in the introduction J/ ip- and tfj'-meson or pair of ^Z^'- mesons can be produced by 
the same mechanism in gluon-gluon interactions. According to dimensional considerations, 
ratios of yields of different meson pairs are determined by the ratios of their wave functions 
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FIG. 3. Angular distribution of different polarization states in gg 2 J ftp reaction at a/I = 10 GeV: 
solid curve — total cross section; dashed curve — both mesons are transversely polarized; dotted curve 
— one meson is transversely polarized, another — longitudinally; dash-dotted curve — both mesons are 
longitudinally polarized. 



at origin and masses: 



m. 



(3) 



The yield of non-identical particles is factor 2 enhanced compared to the yield of identical 
ones. We estimate contribution of the gg — >■ J/t/jt/j' process by substituting {m{J/ip) + 
m{ip'))/2 for m{J/ip) in the matrix element and using the appropriate ip{0) for the second 
quarkonium. For tjj' this value determined from the leptonic width is equal 



V'v,'(r)|,=o = 0-16GeV=^/'. 
Taking numerical values into account one gets approximately 

a{2J/il)) : a{J/il) + : a{2'4j') ~ 1 : 1/2 : 1/12. 



(4) 



(5) 



Noticing that branching fraction oit/j' ^ J/t/j + X decay is about 56% one concludes that 
feed down from this excited state is nearly 30%. 
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Determination of feed down from P-wave states like Xc requires calculation of another 
matrix element with 2 P-wave particles in the final state. Meanwhile production of Xc + 
states is suppressed by C-parity conservation and feed down from Xc pairs is suppressed by 
branching squared. Naive estimation of feed down from the latter process leads to the value 
of about 6% in low invariant mass region. 

As noted above, near the threshold one of the sources of J/^^-pair is gluonic production 
of cc and cc diquarks. This reaction is described by the diagrams similar to those shown in 



figiH but with changed color factors and value of wave function at origin. The 



ast parameter 



can be calculated in the framework of potential model. According to paper 20|| its numerical 
value is equal 

^HWL=o = 0-15GeV3'''- (6) 

In figJUthe cross section of double diquark production in gluonic interaction as a function of 
its invariant mass is shown. Below the threshold of two Scc-baryons production at irigg ^ 7 
GeV we expect the transition of this state into a pair of J/-?/;- mesons or a tetraquark, that 
can be observed as a peak in the J/ipJ/ip mass spectrum. 

III. PRODUCTION AT LHC 

To calculate the cross section of charmonia pairs production in hadronic experiments 
one needs to convolute partonic cross section presented in the previous section with the 
distribution functions of partons in the initial hadrons: 

da {pp ^ 2J/ip + X) = J dxidx2fi{xi) fj {X2) da(ij 2J/ip), (7) 

where xi_2 are the momentum fractions carried by partons. As already mentioned, gluon 
interaction gives main contribution at the LHC energies, so it is sufficient to take only gluon 
fusion into account. Typical values of x are small, xi_2 ~ 10^'^. According to our calculations 
total cross section of the direct J/ ip pairs production in proton-proton interaction at ^/s = 7 
TeV is equal to 

a {pp 2J/4j + X) = 18 nb. 



This value agrees well with previous results ({g], for example). Leading order expression for 
the running strong coupling constant as{rhT) was used to obtain this value. It is clear that 
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since the cross section of the gg 2J/tp process is proportional to final result depends 



22| parametrization at 



strongly on the choice of this constant. For gluon PDFs CTEQ5L 
TTiT scale was used as the most natural one for the leading order QCD process. 

It worth mentioning that simple expression ([7]) does not account for the transverse mo- 
mentum of gluons induced by the radiation in the initial state (ISR). This effect obviously 
does not affect value of the total cross section, but may be significant for the transverse 
momentum distribution of the final particles as well as for the cross sections in different 



detectors. In our work ISR is taken into account by the Pythia 6.4 MC generator 21 1. 
Standard tune of the version 6.4.25 was used. 

To obtain predictions for specific experiment it is necessary take kinematical constraints 
imposed by the detector into account. For instance, at the LHCb experiment main limitation 
is imposed on the rapidity of the final charmonium: 2 < y < 4.5, while there is virtually no 
cutoff on the transverse momentum. Fig. [5] shows rapidity distribution of the J/ip-meson 
at the ^/s = 7 TeV energy. Solid and dashed lines correspond to the transverse motion of 
partons accounted for and neglected, respectively. Cross section of the J/ip pair production 
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FIG. 5. Rapidity distribution of the final charmonium in the reaction 2J/ip+X. Solid and dashed lines 
show the results obtained with the initial state radiation accounted for and in the collinear approximation, 
respectively. 

with the 2 < y < A.5 restriction imposed is equal 

<HCb (pp^2J/^ + X) = 3.2nb 
if the initial state radiation is neglected, and 

cTLHCb {pp 2J/^ + X) = 3.1 nb 

if accounted for. One can see that refusal from the collinear approximation affects cross 
section of the double J/ ip production in LHCb very slightly. This guaranties that uncertainty 
induced by the ISR effect on the fraction of total cross section selected by the cutoff is small 
compared to uncertainty in value of the total cross section itself. The latter dues to a^, 
i^j/ipiS^) and rric values and reaches 30%. 

The case is completely different for the detectors examining central region and limiting 
minimal transverse momentum (such as ATLAS and CMS). The restrictions on the trans- 
verse momentum and rapidity of the final charmonia are: 

Pr>5GeV, |y| < 2.5. 

Fig. [6] shows distributions over the transverse momentum of one of the J/t/)- mesons obtained 
in the collinear approximation and taking the ISR into account. It can be seen that account- 
ing for the ISR leads to a substantial broadening of the transverse momentum distribution 
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FIG. 6. Transverse momentum distribution of the single J/i/i-meson obtained with the initial state radiation 
accounted for (solid line) and in the coUinear approximation (dashed line). 

as compared to the collinear approximation. According to our calculations cross section 
of the J/ip pair production in the collinear approximation taking the ATLAS cutoffs into 
account equals 



while accounting for the transverse motion of gluons lead (with the same cutoffs) to the 
value 



Opposite to the situation in LHCb, this value is very mo del- dependent. Moderate ISR de- 
creases fraction of events selected by the cutoff, as for two particles with transverse momenta 
slightly above the cutoff it often moves one of them below the cutoff value. Intensive ISR 
can increase this fraction as momenta of initial gluons more often become compatible with 
the cutoff value and more particles exceed the cutoff. Meanwhile behavior of the gluon's 
distribution over in the high p^ region is worst known. Values of the cross section under 
ATLAS conditions can differ in several times depending on the tune of Pythia used. 

To demonstrate importance of the initial state radiation we present the py-distribution of 
the J/V'-pair (fig IT]). In the collinear approximation this distribution is obviously described 
by the 5-function. 

Solid curve in fig. [8] shows distribution over the invariant mass of the J / ip pair. One can 
see that without cutoff in transverse momentum main part of events is concentrated in the 



(^atlaT (pP 2 J/^ + X) = 0.09 nb, 



(TATLAS {PP ^ 2 J/^ + X)= 0.06 nb. 
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FIG. 7. Distribution over the transverse momentum of the J /ip pair. 
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FIG. 8. Distribution over the invariant mass of charmonium pairs in the pp — )■ J/^ + X (sohd curve) and 
pp — >■ J/ip + tp' (dashed curve) reaction. 

low invariant mass region. Taking the transverse motion of the gluons into account virtually 
does not affect shape of the invariant mass distribution, as expected. Dashed curve in fig. [8] 
corresponds to the production of the J/ip + final state. This distribution obviously starts 
at bigger invariant mass value than those of the J/ip-paiT. Meanwhile invariant masses of 
J/ip-paiis in which one of J/^/'-mesons originate from the ip' decay {ip' — )■ J/ipTcir mostly) are 
distributed similarly to direct J/ ip-paiis. Both distributions and their total are presented in 
fig. [9l Kinematical cutoff corresponding the LHCb detector (2 < yj/^ < 4.5) is applied to 
events saturating this distribution to compare with the first LHCb data jl5| . 

Cutoff 2 < yj/^ < 4.5 applied to the pp J/ipip' + X, ip' ^ J/ipim events simulated 
in Pythia generator selects approximately 16% of the total feed down cross section, which 
equals approximately 5.7 nb. This corresponds to the feed down in LHCb equal to 



Jp'—^J/lpTTTT 

LHCb 



{pp J/ipip' + X) = 1.0 nb. 
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FIG. 9. Distribution over the invariant mass of J/ij} pairs produced directly (dashed curve), with one J/il) 
originated from -0' decay (dotted curve) and their total (solid curve). LHCb kinematical cutoff is applied. 
Points with error bars correspond to the first LHCb data (isj . 



Together with the direct cross section this leads to the value 



^LHCb {PP ^ 2 J/^ + X) = 4.1 ± 1.2 nb. 



15| is 5.6 ± l.lnb, which is in a good agreement with the 



Experimental value reported in 
prediction. 

We make predictions for different kinematical distributions which will be measured by 
LHCb later at bigger statistics. Figures fTOl and fTTl show distributions over transverse momen- 
tum of a single J /ip from a pair and over transverse momentum of the whole pair respectively. 
Contributions from both direct production and from ip' decays are shown. One sees that 
shape of both contributions are similar. Fig. [12] shows distribution over the rapidity of one 
of the J /if) mesons from a pair. Smooth left border of this distribution dues to the fact that 
two J/ip mesons in a pair are generally close in rapidity. So when one of them fails cutoff 
in rapidity, second one is rejected too, even if it is in the allowed range. 



IV. TETRA-C-QUARK 

The most interesting phenomenon in the region of small invariant masses of J/^/'-pairs 
is formation of a tetraquark built of 2 c-quarks and 2 c-quarks and its decay into a J/ip- 
pair. Similar states, for example exotic mesons X(3872), X(3940), y(3940) and y(4140), 
have already been observed in i?-meson decays and J/?/'^, J/ipu invariant mass spectrums 
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FIG. 10. Distribution over the transverse momentum of a single J/f/i-meson from the pair produced directly 
(dashed curve), with one J/^ originated from i/j' decay (dotted curve) and their total (solid curve). LHCb 
kinematical cutoff is applied. 
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FIG. 11. Distribution over the transverse momentum of a J/?/;- meson pair produced directly (dashed curve), 
with one J/^ originated from ip' decay (dotted curve) and their total (solid curve). LHCb kinematical cutoff 
is applied. 



Id, 



11| . These states have small widths and can not be described in usual quark-antiquark 



scheme. It turns out, however, that all these mesons can be described in terms of tetraquarks 
[c^ls [c^ls, where q is u-, d- or s-quark . In the following we consider the q = c case, 
that is a tetraquark T^c built from four charm quarks — [cc] [cc] . 

To estimate parameters of these states we use diquark model of tetraquark that is we 
assume that it is built from two almost point-like diquarks which interact with each other. 
In this case Pauli exclusion principle imposes strong restrictions on the quantum numbers 
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FIG. 12. Distribution over the rapidity of a single J/^/;-meson of the pair produced directly (dashed curve), 
with one J/tp originated from ip' decay (dotted curve) and their total (solid curve). LHCb kinematical cutoff 
is applied. 

of these diquarks. Indeed, the angular momentum of quarks in the ground state of the 
diquark should be equal to zero. In order to attract to each other and build a tetraquark 
state both diquark and antidiquark should be in the triplet color state. From the Pauli 
exclusion principle it immediately follows that the total spin of both diquarks must be equal 
to unity. Indeed coordinate component of wave function of these diquarks is symmetric, color 
component — antisymmetric, so spin component of wave function has to be symmetric. Mass 
of a diquark with such quantum numbers can be determined by solving a non-relativistic 
Schrodinger equation [20, 23]. Interaction potential is chosen proportional to the potential 
between quark and antiquark in usual charmonium states but with an additional factor 1/2 
which accounts for the difference in color structures. Using mass of the c-quark equal to 
rric = 1.468 GeV and taking into account hyperfine splitting 



Am 



9m? 



[SciSc2)|^h(0)| ~6.4MeV, 



where Sci,2 are quark spin operators and 'ip[cc]{^) is the diquark wave function at origin, 
which value was given in expression ([6]), the following values for diquark mass and mean 
radius can be obtained: 



m[cc] = 3.13 GeV, 



:]) = 0.523 fm. 



The following consideration assumes that both diquarks are point-like. In this case 



15 



tetraquark mass can be obtained by solving Schrodinger equation with the potential used in 
meson spectrum calculation. The wave function of tetraquark at origin appears to be equal 

^T,Xr)l=, = QA7GeY''\ (8) 

Without spin-spin interaction mass and mean radius of tetraquark are equal 

Mt4, = 6.12 GeV, (rr^,) = 0.29 fm. (9) 

It is seen that the mean radius of tetraquark and its constituent diquarks are comparable, 
and therefore the values obtained for the masses should be regarded as a rough estimation. 
Finite size of the diquarks can be accounted for with the form factor 

f 

F (r^) = exp <^ 

Our calculations however show that this correction increases diquark mass only slightly, 
approximately by 40 MeV. Due to the spin-spin interaction, described by the potential 



Vss{r) = ^^(SiS,)5(r). 

9 mf 



[cc] 

the state with mass ([9]) splits into scalar, axial and tensor tetraquarks with masses 



J = 


: 


^T4c(0+ + ) 


= 5.97 GeV, 


J = 


1 : 




= 6.05 GeV, 


J = 


2 : 




= 6.22 GeV. 



It is evident that only tensor meson from this list is above the threshold of J/ip-meson 
pair production. Corrections caused by the finite size of diquarks can, in principle, slightly 
increase masses and move the axial meson above this threshold, but its production in two 
gluon channel is forbidden by the charge parity conservation. As for the scalar tetraquark 
its decay into two real J/ ip mesons is kinematically forbidden. One of the final J / t/^- mesons, 
however, could be a virtual particle decaying into |U+/i~-pair. For this reason it is very 
interesting to search for the resonance in J/?/'/i"*"/i~ channel. 

Duality relations can be used to estimate cross section of tetraquark production. Let us 
consider reaction gg — > [ccjg + [ccjs which cross section as a function of invariant mass of 
the system is plotted in figJH Below the threshold of 2 doubly-heavy baryons Sec diquarks 
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FIG. 13. Tetraquark contribution to the distribution over the invariant mass of the J/i/'-pair at 7TeV centre 
of mass energy. 

can form a tetraquark T^c consisting of 4 c-quarks with subsequent decay into a. J/ip meson 
pair. So, the following relation should hold: 

S= j dmgg(T(gg ^ X^c^^J/ip) 

= K- J dmgg&igg ^ [ccjg + [cc]^) ^ K ■ 6Aph ■ GeV, (10) 

2Mj/^ 

where K < 1. 

In the current work we use value K = 0.1. Width and height of the peak caused by the 
tetraquark contribution is determined by its own width A ~ 0.1 MeV. However it is con- 
siderably smaller than experimental resolution of the detector Aexp. ~ 50 MeV. Histogram 
in Fig. [13] shows the distribution over the invariant mass of the J/ V^-pair with expected 
contribution from the tetraquark using experimental resolution for the bin width. 

V. CONCLUSION 

In our paper we use perturbative QCD to calculate the fourth order contributions to 
the J/ip pairs production at LHC According to our estimates the ratio of the cross section 
of this reaction to the cross section of single J / ip-meson production in conditions of the 
LHCb experiment {2 < y < 4.5) is about 4 x 10"^. It should be noted, that comparison 
of the first experimental results at low energies with theoretical predictions shown excess of 
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experimental data by a factor 2 — 3 jst . Cross section predicted in our study is in a much 
better agreement with the recent LHCb result. Meanwhile all contributions accounted for 
in our work lead to kinematical peak corresponding to the first bin in experimental data 
presented. Possible contributions from the diquark-antidiquark production also lay in this 
invariant mass interval. If further experimental analysis eliminates lack of signal in the first 
bin preserving shape of the distribution corresponding increase of the cross section measured 
can amount to 50% — 100%. 

It is well known, that production and subsequent radiative decays of P-wave states Xc 
can give additional contribution to J / meson cross section. Our estimates in the region of 
low invariant masses of the J /ip pairs do not support this view, though it requires additional 
quantitative validation. Contribution form Xc states can however become significant at large 
invariant masses and it is important to distinguish it from the color-octet one. The latter 
also is hardly noticeable in the kinematical region considered in our article. Other possi- 
ble source of enhancement is inclusive production of J/^^-pairs accompanied by additional 
gluons. Despite the fact that this process is described in the higher orders of perturbation 
theory, the number of color degrees of freedom in the final state increases greatly. It does 
not seem to be possible to distinguish experimentally exclusive J/^/'-pair production and 
reactions with emission of additional gluons. One can only estimate possible enhancement 
caused by contributions from higher orders of QCD. Thus accurate measurement of the cross 
section of double J/ip production is needed. 

Another interesting question considered is the possibility of observation of a new exotic 
meson built from 2 c-quarks and 2 c-quarks. Our estimates show that if diquarks [cc] and 
[cc] are both in color-triplet states than in addition to attraction of quarks in diquark mutual 
attraction of diquarks exist. Rough estimation of masses of resulting tetraquarks is presented 
in our paper and in is shown that at least one of these states (tensor tetraquark) can be 
observed in J/ipJ/tp or J/ipfi'^fi~ modes. 
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